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ABSTRACT 

We present the results of extensive mult i- waveband monitoring of the blazar 
3C 279 between 1996 and 2007 at X-ray energies (2-10 keV), optical R band, 
and 14.5 GHz, as well as imaging with the Very Long Baseline Array (VLBA) 
at 43 GHz. In all bands the power spectral density corresponds to "red noise" 
that can be fit by a single power law over the sampled time scales. Variations in 
flux at all three wavebands are significantly correlated. The time delay between 
high and low frequency bands changes substantially on time scales of years. A 
major multi-frequency flare in 2001 coincided with a swing of the jet toward a 
more southerly direction, and in general the X-ray flux is modulated by changes 
in the position angle of the jet near the core. The flux density in the core at 43 
GHz — increases in which indicate the appearance of new superluminal knots — is 
significantly correlated with the X-ray flux. 
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We decompose the X-ray and optical light curves into individual flares, finding 
that X-ray leads optical variations (XO) in 6 fiares, the reverse occurs in 3 fiares 
(OX), and there is essentially zero lag in 4 fiares. Upon comparing theoretical 
expectations with the data, we conclude that (1) XO fiares can be explained 
by gradual acceleration of radiating electrons to the highest energies; (2) OX 
fiares can result from either light-travel delays of the seed photons (synchrotron 
self-Compton scattering) or gradients in maximum electron energy behind shock 
fronts; and (3) events with similar X-ray and optical radiative energy output 
originate well upstream of the 43 GHz core, while those in which the optical 
radiative output dominates occur at or downstream of the core. 

Subject headings: galaxies: quasars: general — galaxies: quasars: individual 
(3C 279) — physical data and processes — X-rays: observations — radio contin- 
uum: galaxies 



Introduction 



Blazars form a subclass of active galactic nuclei (AGN) characterized by violent variabil- 
ity on time scales from hours to years across the electromagnetic spectrum. It is commonly 
thought that at radio, infrared, and optical frequencies the variable emission of blazars 



originates in r elativistic jets and is synchrotron in nature ( llmpey &: Neugenbauerl Il988 



Marscheiill998l ). The X-ray emission is consistent with inverse Compton (IC) scattering of 
these synchrotron photons, although seed photons from outside the jet cannot be excluded 



(iMause et al 



2000 



1996 



Sikora et al. 



Ronianova &: Lovelace 19971: Coppi fc Aharonian 1999 : Blazejowski et al. 



2001 



Chiang fc Bottcher 



2002 



Arbeiter et al 



20051 ). The models may 



be distinguished by measuring time lags between the seed-photon and Compton-scattered 
fiux variations. Comparison of the amplitudes and times of peak fiux of fiares at different 
wavebands is another important diagnostic. For this reason, long-term multi-frequency mon- 
itoring programs are crucially important for establishing a detailed model of blazar activity 
and for constraining the physics of relativistic jets. Here we report on such a program that 
has followed the variations in emission of the blazar 3C 279 with closely-spaced observations 
over a time span of ~10 years. 



The quasar 3C 279 (redshift=0.538; iBurbidge fc Rosenburd 1 19651 ) is one of the most 
prominent blazars owing to its high optical polarization and variability of fiux across the 
electromagnetic spectrum. Very long baseline interferometry (VLBI) reveals a one-sided 
radio jet featuring bright knots (comp onents) that are "ejected" from a bright, presumably 
stationary "core" Jjorstad et al.ll2005h . The measured apparent speeds of the knots observed 
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in the past range from 4c to 16c (jJorstad et al.ll2004l ). superluminal motion that results from 
relativistic bulk velocities and a small angle between the jet axis and line of sight. Relativistic 
Doppler boosting of the radiation increases the apparent luminosity to ~ 10^ times the value 
in the rest frame of the emitting plasma. 

Characterization of the light curves of 3C 279 includes the power spectral density (PSD) 
of the variability at all different wavebands. The PSD corresponds to t he power in the vari 



ability of emission as a function of time scale. iLawrence et al.l (119871 ) and lMcHardy fc Czerny 



(119871 ) have found that the PSDs of many Seyfert galaxies, in which the continuum emission 
comes mainly from the central engine, are simple "red noise" power laws, with slopes be- 
tween — 1 and —2. More recent studies indicate that some Seyfert galaxies have X-ray PSDs 



(Uttlev et al. 


2002 


Pounds et al. 


2001) 



McHardv et al.l |2004 iMarkowitz et al.l l2003l : lEdelson fc Nandral 11999 



200ll ). This property of Seyferts is similar to that of Galactic black hole X-ray 



binaries (BHXRB) , whose X-ray PSDs contain one or more breaks (IBelloni fc Hasingerlll990 
Nowak et al.lll999l ). However, in blazars most of the X-rays are likely produced in the jets 
rather than near the central engine as in BHXRBs and Seyferts. It is unclear a priori what 
the shape of the PSD of nonthermal emission from the jet should be, a question that we 
answer with the dataset presented here. 

The raw PSD calculated from a light curve combines two aspects of the dataset: (1) the 
intrinsic variation of the object and (2) the effects of the temporal sampling pattern of the 
observations. In order to remove the latter, we apply a Monte-Carlo type alg orithm based 
on the "Power Spectrum Response Method" (PSRESP) of lUttley et al.l (120021 ) to determine 
the intrinsic PSD (and its associated uncertainties) of the light curve of 3C 279 at each of 
three wavebands. Similar complications affect the determination of correlations and time 
lags of variable emission at different wavebands. Uneven sampling, as invariably occurs, 
can cause the correlation coefficients to be artificially low. In addition, the time lags can 
vary across the years owing to physical changes in the source. In light of these issues, we 
use simulated light curves, based on the underlying PSD, to estimate the significance of the 
derived correlation coefficients. 

In §2 we present the observations and data reduction procedures, while in §3 we describe 
the power spectral analysis and its results. In §4 we cross-correlate the light curves to deter- 
mine the relationship of the emission at different wavebands. Finally, in §5 we discuss and 
interpret the results, focusing on the implications regarding the location of the nonthermal 
radiation at different frequencies, as well as the physical processes in the jet that govern the 
development of flux outbursts in blazars. 
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2. Observations and Data Analysis 

2.1. Monitoring of the X-ray, Optical, and Radio Flux Density 

Table [T] summarizes the intervals of monitoring at different frequencies for each of the 
three wavebands in our program. We term the entire light curve "monitor data"; shorter 
segments of more intense monitoring are described below. 

The X-ray light curves are based on observations of 3C 279 with the Rossi X-ray Timing 
Explorer (RXTE) from 1996 to 2007. We observed 3C 279 in 1222 separate pointings with the 
RXTE, with a typical spacing of 2-3 days. The exposure time varied, with longer on-source 
times — typically 1-2 ks — after 1999 as the number of fully functional detectors decreased, 
and shorter times at earlier epochs. For each exposure, we used routines from the X-ray 
data analysis software FTOOLS and XSPEC to calculate and subtract an X-ray background 
model from the data and to model the source spectrum from 2.4 to 10 keV as a power law 
with low-energy photoelectric absorption by intervening gas in our Galaxy. For the latter, 
we used a hydrogen column density of 8 x 10^° atoms cm~^. There is a ~ 1-year gap in 
2000 and annual 8-week intervals when the quasar is too close to the Sun's celestial position 
to observe. A small part of the X-ray data is given in Table [2J The whole dataset will be 
available in the electronic version of ApJ. 

In 1996 December and 1997 January we obtained, on average, two measurements per 
day for almost two months. We refer to these observations as the "longlook" data. Between 
2003 November and 2004 September, we obtained 127 measurements over 300 days (the 
"medium" data). Figure [1] presents these three datasets. The X-ray spectral index a^., 
defined by /x oc z/"^, where /x is the X-ray flux density and v is the frequency, has an 
average value of —0.8 with a standard deviation of 0.2 over the ~ 10 years of observation, 
and remained negative throughout. 

We monitored 3C 279 in the optical R band over the same time span as the X-ray 
observations. The majority of the measurements between 1996 and 2002 are from the 0.3 
m telescope of the Foggy Bottom Observatory, Colgate University, in Hamilton, New York. 
Between 2004 and 2007, the data are from the 2 m Liverpool Telescope (LT) at La Palma, 
Canary Islands, Spain, supplemented by observations at the 1.8 m Perkins Telescope of 
Lowell Observatory, Flagstaff, Arizona, 0.4 m telescope of the University of Perugia Obser- 
vatory, Italy, 0.7 m telescope at the Crimean Astrophysical Observatory, Ukraine, the 0.6 
m SMARTS consortium telescope at the Cerro Tololo Inter-American Observatory, Chile, 
and the 0.7 m Meniscus Telescope of Abastumani Astrophysical Observatory in Abastumani, 
Republic of Georgia. We checked the data for consistency using overlapping measurements 
from different telescopes, and applied corrections, if necessary, to adjust to the LT system. 
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We processed the data from the LT, Perkins Telescope, Crimean Astrophysical Observatory, 
and Abastuma ni Astrophysical Obse r vatory in the same manner, using comparison stars 2, 
7, and 9 from I Gonzalez- Perez et al.l (l200ll ) to determine the magnitudes in R band. The 
frequency of optical measurements over the ~ 10-year span presented here is, on average, 
2-3 observations per week. Over a three-month period between 2005 March and June, we 
obtained about 100 data points, i.e., almost one per day ("longlook" data). Another subset 
("medium") contains ~ 100 points over 200 days between 2004 January and July. Figure [2] 
displays these segments along with the entire 10-year light curve. A small part of the optical 
data is given in Table [31 The whole dataset will be available in the electronic version of ApJ. 

We have compiled a 14.5 GHz light curve (Fig. [3]) with data from the 26 m antenna 
of the University of Michigan Radio Astronomy Obse rvatory. Details of the calibration and 
analysis techniqu es are described in lAller et al. The flux scale is set by observations 

of Cassiopeia A (IBaars et al.l 119771 ). The sampling frequency was usually of order once per 
week. An exception is a span of about 190 days between 2005 March and September when 
we obtained 60 measurements, averaging one observation every ~ 3 days ("medium" data). 
A small part of the radio data is given in Table HI The whole dataset will be available in the 
electronic version of ApJ. 



2.2. Ultra-high Resolution Images with the Very Long Baseline Array 



Starting in 2001 May, we observed 3C 279 with the Very Long Baseline Array (VLBA) 
at roughly monthly intervals, with some gaps of 2-4 months. The sequence of images from 
these data (Fig. H] to Fig. ^ provides a dynamic view of the jet at an angular resolution 
~ 0.1 milliarc s econd s (mas). We processed the data in the same manner as described in 
Jorstad et al. (120051 ) . For epoch s frorn 1995 to 2001, we use the imag es and results of 
Lister et all Jl998h . IWehrle et all J200lh . and Ijorstad et all f l200ll . l2005h . We model the 



brightness distribution at each epoch with multiple cir cular Gaussian c omponents using the 
task MODELFIT of the software package DIFMAP JshepherdI ^99?! ) . At each of the 80 
epochs of VLBA observation since 1996, this represents the jet emission downstream of the 
core by a sequence of knots (also referred to as "components" ) , each characterized by its flux 
density, FWHM diameter, and position relative to the core. Figure [10] plots the distance 
vs. epoch for all components brighter than 100 mJy within 2.0 mas of the core. We use the 
position vs. time data to determine the projected direction on the sky of the inner jet, as 
well as the apparent speeds and birth dates of new superluminal knots. 

We define the inner-jet position angle (PA) ^jet with respect to the core as that of 
the brightest component within 0.1-0.3 mas of the core. As seen in Figure [TTl O^^t changes 
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significantly (~ 80°) over the 11 years of VLBA monitoring. Figure [T2l displays a sampling 
of the VLBA images at epochs corresponding to the circled points in the lower panel of 
Figure [TTJ 

We determine t he apparent speed Pa pp of the moving components using the same pro- 
cedure as defined in IJorstad et al.l (120051 ) . The ejection time Tq is the extrapolated epoch 
of coincidence of a moving knot with the position of the (presumed stationary) core in the 
VLBA images. In o rder to obtain the most accurate values of Tq, given that non-ballistic 
motions may occur (IJorstad et al.l 12004 l2005l ). we use only those epochs when a component 
is within 1 mas of the core, inside of which we assume its motion to be ballistic, ^jet, To, and 
Papp between 1996 and 2007 are shown in Table [51 As part of the modeling of the images, 
we have measured the flux density of the unresolved core in all the images, and display the 
resulting light curve in Figure [T31 



3. Power Spectral Analysis and Results 

For all three wavebands, we used the monitor, medium, and longlook data to calculate 
the PSD at the low, intermediate, and high frequencies of variation, respectively. Since we 
do not have any longlook data at 14.5 GHz, we determine the radio PSD up to the highest 
variational frequency that can be achieved with the medium data. 



3.1. Calculation of the PSD of the Observed Light Curve 



We follow lUttley et al.l (120021 ) to calculate the PSD of a discretely sampled light curve 
f{ti) of length N points using the formula 



FN(iy) 



N 



fiU) cos(27rz/ti 



i=l 



+ 



N 



i=l 



(1) 



This is the square of the modulus of the discrete Fourier transform of the (mean subtracted) 
light curve, calculated for evenly spaced frequencies between i/min and z/max, i-e., i^mm, 2z/min, • 
• ^max- Here, h'^i^=l/T (T is the total duration of the light curve, ^at — ti) and h'jaax=N /2T 
equals the Nyquist frequency z^Nyq- We use the following normalization to calculate the final 
PSD: 



2T 



(2) 



where /i is the average flux density over the light curve. 
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We bin the data in time intervals AT ranging from 0.5 to 25 days, as listed in Table O 
averaging all data points within each bin to calculate the flux. For short gaps in the time 
coverage, we fill empty bins through linear interpolation of the adjacent bins in order to 
avoid gaps that would distort the PSD. We account for the effects of longer gaps, such as 
sun- avoidance intervals and the absence of X-ray data in 2000, by inserting in each of the 
simulated light curves the same long gaps as occur in the actual data. 



3.2. PSD Results: Presence/absence of a Break 



We use a variant of PSRESP (lUttley et al.ll2002l ) to determine the intrinsic PSD of each 



light curve. The method is described in the Appendix. PSRESP gives both the best-fit PSD 
model and a "success fraction" that indicates the goodness of fit of the model. 

The PSDs of the blazar 3C 279 at X-ray, optical, and radio frequencies show red noise 
behavior, i.e., there is higher amplitude variability on longer than on shorter timescales. The 
X-ray PSD is best fit with a simple power law of slope —2.3 ± 0.3, for which the success 
fraction is 45%. The slope of the optical PSD is —1.7 ± 0.3 with success fraction 62%, and 
for the radio PSD it is —2.3 ± 0.5 with success fraction 96%. The observed PSDs and their 
best-fit models are shown in Figure [HI The error bars on the slope represent the FWHM 
of the success fraction vs. slope curve (Figure \TE\i . The rejection confidence, equal to one 
minus the success fraction, is much less than 90% in all three cases (55%, 38%, and 4% in 
the X-ray, optical, and radio wavebands, respectively). This implies that a simple power-law 
model provides an acceptable fit to the PSD at all three wavebands. 

We also fit a broken power-law model to the X-ray PSD, setting the low-frequency 
slope at —1.0 and allowing the break frequency and the slope above the break over a wide 
range of paramet ers (10~^ to 10~^ Hz an d —1.0 to —2.5, respectively) while calculating the 



success fractions ( iMcHardy et al.l (120061 )). Although this gives lower success fractions than 
the simple power-law model for the whole paramater space, a break at a frequency < 10~^ 
Hz with a high frequency slope as steep as —2.4 cannot be rejected at the 95% confidence 
level. 



4. Cross-correlation Analysis and Results 

We employ the discrete cross- correlation function fDCCFi lidelson fc Kroli^ll988h method 
to find the correlation between variations at pairs of wavebands. We bin the light curves of 
all three wavebands in 1-day intervals before performing the cross-correlation. In order to 
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determine the significance of tlie correlation, we perform tlie following steps : 

1. Simulation of M (we use M=100) artificial light curves generated with a Monte-Carlo 
algorithm based on the shape and slope of the PSD as determined using PSRESP for both 
wavebands (total of 2M light curves). 

2. Resampling of the light curves with the observed sampling function. 

3. Correlation of random pairs of simulated light curves (one at each waveband). 

4. Identification of the peak in each of the M random correlations. 

5. Comparison of the peak values from step 4 with the peak value of the real correlation 
between the observed light curves. For example, if 10 out of 100 random peak values are 
greater than the maximum of the real correlation, we conclude that there is a 10% chance 
of finding the observed correlation by chance. Therefore, if this percentage is low, then the 
observed correlation is significant even if the correlation coefficient is substantially lower 
than unity. 

As determined by the DCCF (Figure [T6l) . the X-ray variations are correlated with those 
at both optical and radio wavelengths in 3C 279. The peak X-ray vs. optical DCCF is 0.66, 
which corresponds to a 98% significance level. The peak X-ray vs. radio DCCF is relatively 
modest (0.42), with a significance level of 79%. The radio-optical DCCF has a similar peak 
value (0.45) at a 62% significance level. The cross-correlation also indicates that the optical 
variations lead the X-ray by 20 ± 15 and the radio by 260l^°o days, while X-rays lead the 
radio by 2401^qo days. The uncertainties in the time delays are the FWHM of the peaks in 
the correlation function. 



4.1. Variation of X-ray /optical Correlation with Time 

The X-ray and the optical light curves are correlated at a very high significance level. 
However, the uncertainty in the X-ray-optical time delay is comparable to the delay itself. To 
characterize the variation of the X-ray/optical time lag over the years, we divide both light 
curves into overlapping two-year intervals, and repeat the DCCF analysis on each segment. 
The result indicates that the correlation function varies significantly with time (Fig. [T7|) over 
the 11 years of observation. Of special note are the following trends: 

1. During the first four years of our program (96-97, 97-98, 98-99) the X-ray variations lead 
the optical (negative time lag). 

2. There is a short interval of weak correlation in 1999-2000. 

3. In 2000-01, the time delay shifts such that the optical leads the X-ray variations (positive 
time lag). This continues into the next interval (2001-02). 

4. In 2002-03, there is another short interval of weak correlation (not shown in the figure). 
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5. In the next interval (2003-04), the delay shifts again to almost zero. 

6. Over the next 3 years the correlation is relatively weak and the peak is very broad, 
centered at a slightly negative value. 

This change of time lag over the years is the main reason why the peak value of the overall 
DCCF is significantly lower than unity. We discuss the physical cause of the shifts in cross- 
frequency time delay in §5.2. 



4.2. Correlation of X-ray Flux and Position Angle of the Inner Jet 

We find a significant correlation (maximum DCCF=0.6) between the PA of the jet and 
the X-ray flux (see Figure [18]). The changes in the position angle lead those in the X-ray flux 
by 80 ± 150 days. The large uncertainty in the time delay results from the broad, nearly flat 
peak in the DCCF. This implies that the jet direction modulates rather gradual changes in 
the X-ray flux instead of causing specific flares. This is as expected if the main consequence 
of a swing in jet direction is an increase or decrease in the Doppler beaming factor on a time 
scale of one or more years. 



4.3. Comparison of X-ray and Optical Flares 



We follow IValtaoja et al.l (119991 ) by decomposing the X-ray and optical light curves 
into individual flares, each with exponential rise and decay. Our goal is to compare the 
properties of the major long-term flares present in the X-ray and optical lightcurves. Before 
the decomposition, we smooth the light curve using a Gaussian function with a 10-day 
FWHM smoothing time. We proceed by first fitting the highest peak in the smoothed light 
curve to an exponential rise and fall, and then subtracting the flare thus fit from the light 
curve. We do the same to the "reduced" hght curve, i.e., we fit the next highest peak. 
This reduces confusion created by a flare already rising before the decay of the previous 
flare is complete. We fit the entire light curve in this manner with a number of individual 
(sometimes overlapping) flares, leaving a residual flux much lower than the original flux at 
all epochs. We have determined the minimum number of flares required to adequately model 
the lightcurve to be 19 (X-ray) and 20 (optical), i.e., using more than 19 flares to model the 
X-ray light curve does not change the residual flux significantly. 

Figure [19] compares the smoothed light curves with the summed flux (sum of contribu- 
tions from all the model flares at all epochs). We identify 13 X-ray/optical flare pairs in 
which the flux at both wavebands peaks at the same time within ±50 days. Since both light 
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curves are longer than 4200 days and there are only about 20 significant flares during this 
time, it is highly probable that each of these X-ray/optical flare pairs corresponds to the 
same physical event. There are some X-ray and optical flares with no significant counterpart 
at the other waveband. We note that this does not imply complete absence of flaring activity 
at the other wavelength, rather that the corresponding increase of flux was not large enough 
to be detected in our decomposition of the smoothed hght curve. 

We calculate the area under the curve for each flare to represent the total energy output 
of the outburst. In doing so, we multiply the R-band flux density by the central frequency 
(4.7 X 10^^ Hz) to estimate the integrated optical flux. For each of the flares, we determine 
the time of the peak, width (deflned as the mean of the rise and decay times), and area 
under the curve from the best flt model. Table [7] lists the parameters of each flare pair, 
along with the ratio Cxo of X-ray to optical energy output. The time delays of the flare 
pairs can be divided into three different classes: X-ray signiflcantly leading the optical peak 
(XO, 6 out of 13), optical leading the X-ray (OX, 3 out of 13), and nearly coincident (by 
< 10 days, the smoothing length) X-ray and optical maxima (C, 4 out of 13). The number 
of events of each delay classiflcation is consistent with the correlation analysis (Figure fTTll . 
XO flares dominate during the flrst and last segments of our program, but OX flares occur 
in the middle. In both the DCCF and flare analysis, there are some cases just before and 
after the transition in 2001 when variations in the two wavebands are almost coincident (C 
flares) . 

The value of (xo ~ 1 in 5 out of 13 cases; in one flare pair (xo = 1-4. In all the other 
cases it is less than unity by a factor of a few. In all the C flares (xo ~ 1, while in the 3 OX 
cases the ratio -C 1. In the C pairs, the width of the X-ray flare proflle ~ 2 times that of 
the optical, but in the other events the X-ray and optical widths are comparable. 



4.4. Flare-Ejection Correlation 



The co re region on VLBI ima ges becomes brighter as a new superluminal knot passes 
through it (jSavolainen et al.ll2002l ): hence, maxima in the 43 GHz light curve of the core 
indicate the times of ejection of knots. We flnd that the core (Figure [T3|) and X-ray light 
curves are well correlated (correlation coefficient of 0.6), with changes in the X-ray fiux 
leading those in the radio core by 130+^° days (see Figure 1201) . The broad peak in the cross- 
correlation function suggests that the flare-ejection time delay va r ies oy er a rather broad 
range. This result is consistent with the flnding of iLindfors et al.l (120061 ) that high-energy 
flares generally occur during the rising portion of the 37 GHz light curve of 3C 279. 
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Since some flares can be missed owing to overlapping declines and rises of successive 
events, we can determine the times of superluminal ejections more robustly from the VLBA 
data. Table O lists the ejection times and apparent speeds of the knots identified by our pro- 
cedure. We cannot, however, associate an X-ray flare with a particular superluminal ejection 
without further information, since flux peaks and ejection times are disparate quantities. We 
pursue this in a separate paper that uses light curves at five frequencies between 14.5 and 
350 GHz to analyze the relationship between superluminal knots and flares in 3C 279. 



5. Discussion 



5.1. Red Noise Behavior and Absence of a Break in the PSD 



The PSDs at all three wavebands are best fit with a simple power law which corresponds 
to red noise. The red noise nature — greater amplitudes on shorter time scales — of the flux 
variations at all three wavebands revealed by the PSD analysis is also evident from visual 
inspection of the light curves of 3C 279 (Figures [H [21 and [3]). 

The PSD break frequency in BHXRBs and Seyferts scales with the rnass of the black hole 
Juttlev et al.1 12OO2I : iMcHardv etaP booi l20od : iMarkowitz et all I2OO3I : Isdelson fc Nandra 
I999I ). Using the best-fit values and uncertainties in the relation b etwee n break timescale, 
black- hole mass, and accretion rate obtained by iMcHardy et al.l (120061 ). we estimate the 
expected value of the break frequency in the X-ray PSD of 3C 279 to be io~^-^^°-^ Hz, which 
i s just within our d e rived lower limit of 10~ ^'^ Hz. Here we use a black-hole mass of 10^ M© 
( IWoo fc Urryl l2002l : iLiu. Jiang fc Gull2006l) for 3C 279 and a bolom etric luminosity of the 
big blue bump of 4 x 10'^^ ergs/s (IHartman et al.l 119961 ). If we follow iMcHardy et al.l (120061 ) 
and set the low-frequency slope of the X-ray PSD at —1.0 and allow the high-frequency slope 
to be as steep as —2.4, a break at a frequency < 10^^ Hz cannot be rejected at the 95% 
confidence level. An even longer light curve is needed to place more stringent limits on the 
the presence of a break at the expected frequency. 



5.2. Correlation between Light Curves at Different Wavebands 

The cross-frequency time delays uncovered by our DCCF analysis relate to the relative 
locations of the emission regions at the different wavebands, which in turn depend on the 
physics of the jet and the high-energy radiation mechanism(s). If the X-rays are synchrotron 
self-Compton (SSC) in nature, their variations may lag the optical flux change s owing to the 



travel time of the seed photons before they are up-scattered. As discussed in ISokolov et al. 
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( 120041 ). this is an important effect provided tliat tlie angle (6'obs) between tlie jet axis and 
tlie line of sight in the observer's frame is sufficiently small, < 1.2° ± 0.2°, in the case 
of 3C 279, w h ere w e have adopted the bulk Lore n tz fac tor (T = 15.5 ± 2.5) obtained by 
Jorstad et al.l (120051 ) . According to ISokolov et al.l (120041 ) , if the emission region is thicker 
(thickness ~ radius), the allowed angle increases to 2° ± 0.4°. X-rays produced by inverse 
Compton scattering of seed photons from outside the jet (external Compton, or EC, process) 
may lag the low frequency emission for any value of 6'obs between 0° and 90°. However, in 
this case we expect to see a positive X-ray spectral index over a significant portion of a 
f lare if ^ohs is small, and th e flares should be asymmetric, with much slower decay than rise 
(jSokolov &: Marscherll2005l ). This is because the electrons that up-scatter external photons 
(radiation from a dusty torus, broad emission-line clouds, or accretion disk) to X-rays have 
relatively low energies, and therefore have long radiative cooling times. Expansion cooling 
quenches such flares quite slowly, since the EC flux depends on the total number of radiating 
electrons (rather than on the number density), which is relatively weakly dependent on the 
size of the emitting region. 

Time delays may also be produced by frequency stratification in the jet. This occurs 
when the electrons are energized along a surface (e.g., a shock front) and t hen move away at 
a spe ed close to c as they lose energy via synchrotron and IC processes (IMarscher &: Gear 
19851 ). This causes the optical emission to be radiated from the region immediately behind 
the surface, with the IR emission arising from a somewhat thicker region and the radio 
from an even more extended volume. An optical to radio synchrotron flare then begins 
simultaneously (if opacity effects are negligible), but the higher-frequency peaks occur earlier. 
On the other hand, SSC (and EC) X-rays are produced by electrons having a range of 
energies that are most ly lower than those required to produce optical synchrotron emission 
(IMcHardy et al.lll999l ). Hence, X-rays are produced in a larger region than is the case for 
the optical emission, so that optical flares are quenched faster and peak earlier. Flatter PSD 
in the optical waveband is consistent with this picture. 

In each of the above cases, the optical variations lead those at X-ray energies. But 
in majority of the obse rved flares, the revers e is true. This may be explained by another 
scenario, mentioned by iBottcher et al.l (120071 ). in which the acceleration time scale of the 
highest-energy electrons is significantly longer than that of the lower-energy electrons, and 
also longer than the travel time of the seed photons and/or time lags due to frequency 
stratification. In this case. X-ray flares can start earlier than the corresponding optical 
events. 



We thus have a working hypothesis that XO (X-ray leading) flares are governed by 
gradual particle acceleration. OX events can result from either (1) light-travel delays, since 
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the value of ^obs determined by IJorstad et al.l (120051 ) (2.1° ± 1.1°) is close to the required 



range and could have been smaller in 2001-03 when OX flares were prevalent, or (2) frequency 
stratification. One way to test this further would be to add light curves at 7-ray energies, 
as will be possible with the upcoming Gamma Ray Large Area Space Telescope (GLAST). 
If the X and 7 rays are produced by the same mechanism and the X-ray/optical time lag is 
due to light travel time, we expect the 7-ray/optical time lag to be similar. If, on the other 
hand, the latter delay is caused by frequency stratification, then it will be shorter, since IC 7 
rays and optical synchrotron radiation are produced by electrons of similarly high energies. 

If the synchrotron flare and the resultant SSC flare are produced by a temporary increase 
in the Doppler factor of the jet (due to a change in the direction, Lorentz factor, or both), 
then the variations in flux should be simultaneous at all optically thin wavebands. It is 
possible for the C flare pairs to be produced in this way. Alternatively, the C events could 
occur in locations where the size or geometry are such that the time delays from light travel 
and frequency stratification are very short compared to the durations of the flares. 

As is discussed in §4.2, the X-ray/optical time delay varies significantly over the observed 
period (see Table [7] and Fig. [T7j). As we discuss above, XO flares can be explained by 
gradual acceleration of electrons. The switch in the time delay from XO to OX at the onset 
of the highest amplitude multi- waveband outburst, which occurred in 2001 (MJD 2000 to 
2200), might then have resulted from a change in the jet that shortened the acceleration 
time significantly. This could have caused the time delay from light travel time of the 
SSC seed photons and/or the effects of gradients in maximum electron energy to become 
a more significant factor than the acceleration time in limiting the speed at which flares 
developed during the outburst period. However, the flux from the 43 GHz core also reached 
its maximum value in early 2002, and the apparent speed of the jet decreased from ~ 17c 
in 2000 to ~ 4 — 7c in 2001-2003. This coincided with the onset of a swing toward a more 
southerly direction of the trajectories of new superluminal radio knots. We hypothesize 
that the change in direction also reduced the angle between the jet and the line of sight, so 
that the Doppler factor 5 of the jet increased significantly, causing the elevated flux levels 
and setting up the conditions for major flares to occur at all wavebands. The pronounced 
variations in flux during the 2001-2002 outburst cannot be explained solely by fluctuations 
in (5, however, since the time delay switched to OX rather than C. Instead, a longer-term 
switch to a smaller viewing angle would have allowed the SSC light-travel delay to become 
important, causing the switch from XO to OX flares. 
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5.3. Quantitative Comparison of X-ray and Optical Flares 



The relative amplitude of synchrotron and IC flares depends on which physical param- 
eters of the jet control the flares. The synchrotron flux is determined by the magnetic field 
B, the total number of emitting electrons Ng, and the Doppler factor of the flow 6. The IC 
emission depends on the density of seed photons, number of electrons available for scattering 
Ne, and 6. An increase solely in would enhance the synchrotron and EC flux by the same 
factor, while the corresponding SSC flare would have a higher relative amplitude owing to 
the increase in both the density of seed photons and number of scattering electrons. If the 
synchrotron flare were due solely to an increase in B, the SSC flare would have a relative 
amplitude similar to that of the synchrotron flare, since B affects the SSC output only by 
increasing the density of seed photons. In this case, there would be no EC flare at all. Fi- 
nally, if the synchrotron flare were caused solely by an increase in 6, the synchrotron and 
SSC flux would rise by a similar factor, while the EC flare would be more pronounced, since 
the density and energy of the incoming photons in the plasma frame would both increase by 
a factor 6. Table [S] summarizes these considerations. 

The location of the emission region should also have an effect on the multi-waveband 
nature of the flares. The magnetic field and electron energy density parameter Nq both de- 
crease with distance r from the base of the jet: B r^^ and Ai'o ~ r~"; we adopt n = 2 and 
6 = 1 and assume a conical geometry. The cross-sectional radius R of the jet expands with 
r, R (X r. We have performed a theoretical calculation of the energy output of flares that 
includes the dependence on the location of the emission region. We use a computer code 
that calculates the synchrotron and SSC radiation from a source with a power-law energy 
distribution of electrons A^(7) = A^o7~* within a range 7min to 7max5 where 7 is the energy 
in units of the electron rest mass. We introduce time variability of the radiation with an 
exponential rise and decay in B and/or No with time. In addition, we increase 'jmax with 
time in some computations in order to simulate the gradual acceleration of the electrons. 
The synchrotron emission coefficient is given by 



consistent with the mean X-ray spectral index. The critical frequency, near which most 
of the synchrotron luminosity occurs, is given by Uc = fc'7^, while the synchrotron energy 
loss rate is given by dE/dt = —k'-f'^. Both k and k' are functions of B and are given by 
k = 1.3 X 10~^-B^ and k' = 4.2 x 10^5. The inverse-Compton (SSC in this case) emission 
coefficient is given by 




(3) 



where Ks is the modified Bessel function of the second kind of order |. We adopt s = 2.5 



3 
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Ju = [ [ '^Ml^^)Rcrie,,ef,^)Ni^)dJdui, (4) 

where we approximate that the emission/scattering region is uniform (i.e., we ignore fre- 
quency stratification) and spherical with radius R. The Compton cross-section a is a function 
of 7 as well as the incident (ui) and scattered (uf) frequencies of the photons: 

a{ui,Uf,-f) = A^r--L[8 + 2x-x2 + 4xln(|)], (5) 

where x = u f / [ui'-f'^) and cr-r is the Thompson cross-section. 

In Figure [211 the top left panel shows the synchrotron and SSC flares at a given distance 
r from the base of the jet, where r is a constant times a factor af^c (ctfac=l)- The right panel 
shows the same at a farther distance (ofac = 5). The flares are created by an exponential 
rise and decay in the B field, while the value of 7max is held constant. The total time of the 
flares is fixed at 10"^ seconds (120 days). We can see from Figure EH that at larger distances 
(right panel) the SSC flare has a lower amplitude than the synchrotron flare even though 
at smaller values of r (left panel) they are comparable. The bottom panels of Figure [21] 
are the same as above except that the value of 'jmax increases with time, as inferred from 
the forward time delay that occurs in some events. We see the same effect as in the top 
panels. In Figure [221 the top and bottom panels show similar results, but in these cases the 
flares are created by an exponential rise and fall in the value of A^o- Again, the results are 
similar. Figure [231 shows segments of the actual X-ray and optical light curves, which are 
qualitatively similar to the simulated ones. The energy output of both synchrotron and SSC 
flares decreases with increasing r, but more rapidly in the latter. As a result, the ratio of 
SSC to synchrotron energy output (xo decreases with r. The case (xo <^ 1 is therefore a 
natural consequence of gradients in B and Nq. 

From Table [71 we can see that in 7 out of 13 of the flare pairs (xo <C 1. In one pair, 
(xo > 1; and in all other pairs, (xo ~ 1- Our theoretical calculation suggests that the 
pairs where the ratio is less than 1 are produced at a larger distance from the base of the 
jet than those where the ratio > 1. The size of the emission region should be related to the 
cross-frequency time delay, since for all three explanations of the lag a larger physical size of 
the emission region should lead to a longer delay. We can then predict that the X-ray/optical 
time delay of the latter flares should be smaller than for the pairs with ratio < 1. Indeed, 
inspection of Table [71 shows that, for most of the pairs, shorter time delays correspond to 
larger (xo, as expected. The smaller relative width of the optical C flares supports the 
conclusion that these occur closer to the base of the jet than the other flare pairs. 
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6. Conclusions 

This paper presents well-sampled, decade-long light curves of 3C 279 between 1996 
and 2007 at X-ray, optical, and radio wavebands, as well as monthly images o btained with 



the V LB A at 43 GHz. We have applied an algorithm based on a method by lUttley et al. 



(I2OO2I ) to obtain the broadband PSD of nonthermal radiation from the jet of 3C 279. Cross- 
correlation of the light curves allows us to infer the relationship of the emission across different 
wavebands, and we have determined the significance of the correlations with simulated light 
curves based on the PSDs. Analysis of the VLBA data yields the times of superluminal 
ejections and reveals time variations in the position angle of the jet near the core. We have 
identified 13 associated pairs of X-ray and optical flares by decomposing the light curves into 
individual flares. Comparison of the observed radiative energy output of contemporaneous X- 
ray and optical flares with theoretical expectations has provided a quantitative evaluation of 
synchrotron and SSC models. We have discussed the results by focusing on the implications 
regarding the location of the nonthermal radiation at different frequencies, physical processes 
in the jet, and the development of disturbances that cause outbursts of flux density in blazars. 
Our main conclusions are as follows: 

(1) The X-ray, optical, and radio PSDs of 3C 279 are of red noise nature, i.e., there is higher 
amplitude variability at longer time scales than at shorter time scales. The PSDs can be 
described as power laws with no significant break, although a break in the X-ray PSD at a 
variational frequency < 10"^ Hz cannot be excluded at a 95% confidence level. 

(2) X-ray variations correlate with those at optical and radio wavebands, as expected if 
nearly all of the X-rays are produced in the jet. The X-ray flux correlates with the projected 
jet direction, as expected if Doppler beaming modulates the mean X-ray flux level. 

(3) X-ray flares are associated with superluminal knots, with the times of the latter indicated 
by increases in the flux of the core region in the 43 GHz VLBA images. The correlation has 
a broad peak at a time lag of I30I45 days, with X-ray variations leading. 

(4) Analysis of the X-ray and optical light curves and their interconnection indicates that 
the X-ray flares are produced by SSC scattering and the optical flares by the synchrotron 
process. Cases of X-ray leading the optical peaks can be explained by an increase in the time 
required to accelerate electrons to the high energies needed for optical synchrotron emission. 
Time lags in the opposite sense can result from either light-travel delays of the SSC seed 
photons or gradients in maximum electron energy behind the shock fronts. 

(5) The switch to optical-leading flares during the major multi-frequency outburst of 2001 
coincided with a decrease in the apparent speeds of knots from 16-17c to 4- 7c and a swing 
toward the south of the projected direction of the jet near the core. This behavior, as well 
as the high amplitude of the outburst, can be explained if the redirection of the jet (only a 
l°-2° change is needed) caused it to point closer to the line of sight than was the case before 
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and after the 2001-02 outburst. 

(6) Contemporaneous X-ray and optical flares with similar radiative energy output originate 
closer to the base of the jet, where the cross-section of the jet is smaller, than do flares in 
which the optical energy output dominates. This is supported by the longer time delay in 
the latter case. This effect is caused by the lower electron density and magnetic field and 
larger cross-section of the jet as the distance from the base increases. 

Further progress in our understanding of the physical structures and processes in com- 
pact relativistic jets can be made by increasing the number of wavebands subject to intense 
monitoring. Expansion of such monitoring to a wide range of 7-ray energies will soon be 
possible when GLAST scans the entire sky several times each day. When combined with 
similar data at lower frequencies as well as VLBI imaging, more stringent tests on models 
for the nonthermal emission from blazars will be possible. 
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Astronomy Observatory, a facility of the National Science Foundation operated under coop- 
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Roque de los Muchachos of the Instituto de Astrofisica de Canarias, with financial support 
from the UK Science and Technology Facilities Council. 
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A. Power Spectrum Response Method (PSRESP) 



In this study, we determ ine the shape and slope of the PSD of the light curves of 3C 279 
using the PSRESP method (lUttley et al.l 120021 ). This involves the following steps: 

1. Calculation of the PSD of the observed light curve (PSDobs) with formulas (1) and (2). 

2. Simulation of M artificial light curves of red noise nature with a trial shape (simple power 
law, broken power law, bending power law, etc.) and slope. We use M = 100. 

3. Resampling of the simulated light curves with the observed sampling function. 

4. Calculation of the PSD of each of the resampled simulated light curves (PSDsim.i, i=l, 
M) . The resampling with the observed sampling function (which is irregular) adds the same 
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distortions to the simulated PSDs that are present in the real PSD (PSDobs) 
5. Calculation of two functions similar to 



(AP5Aim) 

and 



Xohs - 2^ (AP.Qn_.__\2 ^ ^ 



2^ (AP.^n.A2 ' y^^) 



Sim J 



where PSDsim is the average of (PSDsim,i) and APSDsim is the standard deviation of (PSDsim,i), 
with i=l, M. 

6. Comparison of Xobs with the Xdist distribution. Let m be the number of Xdisti which 
Xobs is smaller than Xdisti • Then m/M is the success fraction of that trial shape and slope, 
a measure of its success at representing the shape and slope of the intrinsic PSD. 

7. Repetition of the entire procedure (steps 2 to 6) for a set of trial shapes and slopes of 
the initial simulated PSD to determine the shape and slope that gives the highest success 
fraction. We scan a range of trial slopes from —1.0 to —2.5 in steps of 0.1 for the simple 
power-law fit. 

We perform a few additional steps to overcome the distorting effects of finite length and 
discontinuous sampling of the light curves. These steps are implicitly included in the light 
curve simulation (step 2). The light curve of an astronomical source is essentially infinitely 
long, but we have sampled a 10-year long interval of it and are calculating the PSD based 
on that interval. As a result, power from longer (than observed) time scales "leaks" into 
the shorter time scales and hence distorts the observed PSD. This effect, called "red noise 
leak" (RNL), can be accounted for in PSRESP. We overcome this by simulating light curves 
that are more than 100 times longer than the observed light curve. As a consequence, the 
resampled simulated hght curves are a small subset of the originally simulated ones, and 
similar RNL distortions are included in PSDsim.i that are present in PSDobs- 

On the other hand, if a light curve is not continuously sampled, power from frequencies 
higher than the Nyquist frequency (i^Nyq) is shifted or "aliased" to frequencies below i^Nyq- 
The observed PSD in that case will be distorted by the aliased power, which is added to 
the observed light curve from timescales as small as the exposure time {Te^p) of the obser- 
vation (about 1000 seconds for the X-ray light curve). Ideally, we should account for this 
by simulating light curves with a time-resolution as small as 1000 seconds so that the same 
amount of aliasing occurs in the simulated data. T his involves excessive computing time for 



decade-long light curves. To avoid this, we follow lUttley et al.l (120021 ) by simulating light 



curves with a resolution lOTe^^p . To calculate the aliasing power from timescales from Texp 
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to lOTexp, we use an analytic approximation of the level of power added to all frequencies by 
the aliasing, given by 



(27cxp ) 

P{iy)du. (A3) 



^Nyq ^min J u^iyq 

- exp J 



We use PSRESP to account for aliasing at frequencies lower than (lOTe; ^ 



We also add Poisson noise to the simulated light curves: 

P ■ - ^"-^'^^^ (A4) 

noise J, J-/- N ; \-^^) 

where a(i) are observationa l uncertainties. For details of the noise processes, readers should 
refer to luttlev et all J2002h 

The goal of adding the noise and resampling with the observed sampling function is 
to simulate a dataset that has the same properties, including the imperfections, as the 
observed one. This provides a physically meaningful comparison of the observed PSD with 
the distribution of the simulated PSDs. 
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Table 1: Start and end times of observations presented in this study. 





X- 


■ray 


Optical 


Radio 


Dataset 


Start 


End 


Start End 


Start End 


Longlook 


Dec 96 


Jan 97 


Mar 05 Jun 05 




Medium 


Nov 03 


Sep 04 


Jan 04 Jul 04 


Mar 05 Sep 05 


Monitor 


Jan 96 


Jun 07 


Jan 96 Jun 07 


Jan 96 Sep 07 
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Table 2. X-ray (2-10 keV) light curve data from RXTE (first 5 rows). 



MJD^ 


Exp. (s) 


Flux^ 


Error 




Err 


count rate^ 


Err 


105.1275 


688 


1.314C-11 


1.818e-12 


0.466 


0.226 


2.766e+00 


1.722e-01 


106.2820 


720 


1.343c- 11 


1.648c- 12 


0.489 


0.206 


2.810e+00 


1.602e-01 


107.2621 


672 


1.150e-ll 


1.829e-12 


0.425 


0.246 


2.393e+00 


1.625e-01 


108.2144 


672 


1.299e-ll 


1.568e-12 


0.544 


0.221 


2.747e+00 


1.670e-01 


109.2619 


624 


1.173e-ll 


1.440e-12 


0.597 


0.245 


2.474e+00 


1.688e-01 



^MJD = Julian date minus 2450000 

^Units: erg cm~^ 

^ "Energy" spectral index 

^Photon counts s"^ per PCU detector 



Table 3. Optical (R Band) light curve data (first 5 rows). 



MJDi 


mag 


err 


Flux Density^ 


err 


46.9514 


14.383 


0.007 


5.837 


0.038 


71.7330 


15.170 


0.032 


2.827 


0.082 


86.6356 


15.386 


0.032 


2.317 


0.067 


91.6813 


14.740 


0.032 


4.201 


0.122 


100.6226 


14.878 


0.032 


3.700 


0.107 



^MJD = Juhan date minus 2450000 
^ Units: mJy 
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Table 4. Radio (14.5 GHz) light curve data (first 5 rows). 



MJD^ 


Flux Density (Jy) 


error 


97.00 


17.51 


0.20 


98.00 


17.74 


0.19 


101.00 


18.62 


0.15 


103.00 


17.54 


0.44 


104.00 


18.28 


0.22 



^MJD = Julian date minus 
2450000 
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Table 5: Ejection times, apparent speeds, and position angle of superluminal knots. 



Knot 




To 






To (MJD^ ) 




9 (deg)2 


C8 


1996. 


,09± 


0. 


,10 


115± 36 


5.4± 0.7 


-130± 3 


C9 


1996. 


.89± 


0. 


.12 


407± 44 


12.9± 0.3 


-131± 5 


CIO 


1997. 


.24± 


0. 


.16 


536± 58 


9.9± 0.5 


-132± 6 


Cll 


1997. 


.59± 


0. 


.11 


662± 40 


lO.li 1.2 


-135± 4 


C12 


1998. 


,56± 


0. 


,09 


1016d= 33 


16.9± 0.4 


-129± 3 


C13 


1998. 


,98± 


0. 


,07 


1174± 26 


16.4± 0.5 


-130± 4 


C14 


1999. 


,50± 


0. 


,09 


1360± 33 


18.2± 0.7 


-135± 6 


C15 


1999. 


,85± 


0. 


,05 


1487± 18 


17.2± 2.3 


-131± 7 


C16 


2000. 


,27± 


0. 


,05 


1642± 18 


16.9± 3.5 


-140± 8 


C17 


2000. 


.96± 


0. 


.12 


1895± 44 


6.2± 0.5 


-133± 12 


C18 


2001. 


.40± 


0. 


.16 


2054± 58 


4.4± 0.7 


-150± 8 


C19 


2002. 


.97± 


0. 


.12 


2648± 44 


6.6± 0.6 


-133± 7 


C20 


2003. 


.39± 


0. 


.10 


2781± 44 


6.0± 0.5 


-155± 10 


C21 


2004. 


,75± 


0. 


,05 


3280± 18 


16. 7± 0.3 


-147± 7 


C22 


2005. 


,18± 


0. 


,06 


3434± 22 


12.4± 1.2 


-102±17 


C23 


2006. 


,41± 


0. 


,15 


3888± 55 


16.5± 2.3 


-114± 5 



iMJD = Julian date minus 2450000 

^Average position angle of knot within 1 mas of the core. 
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Table 6: Parameters of the light curves for calculation of PSD. 





Dataset 


T (days) 


AT (days) 


log(/mm) 


log(/max) 


■^points 




Longlook 


55.0 


0.5 


-6.67 


-4.93 


111 


X-ray 


Medium 


301.0 


5.0 


-7.40 


-5.93 


127 




Monitor 


4150.0 


25.0 


-8.55 


-6.63 


1213 




Longlook 


86.0 


1.0 


-6.86 


-5.23 


94 


Optical 


Medium 


185.0 


5.0 


-7.17 


-5.92 


77 




Monitor 


4225.0 


25.0 


-8.55 


-6.63 


995 


Radio 


Medium 


189.0 


4.0 


-7.18 


-5.83 


59 




Monitor 


3984.0 


25.0 


-8.53 


-6.63 


609 
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Table 7: Total energy output (area) and widths of flare pairs. 



ID 


X- 


-ray 




Optical 




AT 


TDC 






Time^ 


Area ^ 


Width^ 


Timc^ 


Area^ 


Width^ 


(days) 








1 


119 (1996.10) 


33.2 


6.0 


134 (1996.14) 


76.3 


12.0 


-15 


xo 





44 


2 


717 (1997.74) 


371.8 


90.0 


744 (1997.81) 


929.5 


97.5 


-27 


xo 





40 


3 


920 (1998.30) 


172.8 


50.0 


944 (1998.36) 


278.2 


35.0 


-24 


xo 





62 


4 


1050 (1998.65) 


338.7 


70.0 


1099 (1998.79) 


502.8 


57.5 


-49 


xo 





67 


5 


1263 (1999.24) 


160.7 


77.5 


1266 (1999.24) 


168.9 


42.5 


-3 


c 





95 


6 


1509 (1999.91) 


699.8 


202.5 


1517 (1999.93) 


715.4 


90.0 


-8 


c 





98 


7 


2045 (2001.38) 


103.7 


60.0 


2029 (2001.33) 


568.3 


65.0 


16 


ox 





18 


8 


2151 (2001.67) 


453.6 


62.5 


2126 (2001.60) 


1188.3 


57.5 


25 


ox 





38 


9 


2419 (2002.40) 


193.5 


80.0 


2422 (2002.41) 


222.6 


40.0 


-3 


c 





87 


10 


3185 (2004.50) 


191.8 


92.5 


3191 (2004.52) 


218.6 


55.0 


-6 


c 





88 


11 


3416 (2005.13) 


217.7 


70.0 


3444 (2005.21) 


198.7 


50.0 


-28 


xo 


1 


09 


12 


3792 (2006.16) 


525.3 


120.0 


3814 (2006.22) 


375.6 


67.5 


-22 


xo 


1 


40 


13 


4035 (2006.83) 


324.2 


105.0 


4008 (2006.76) 


1068.4 


75.0 


27 


ox 





30 



1 Units: MJD (yr) 
^Units: 10~^ erg cm~^ 
^ Units: days 

^Ratio of X-ray to optical energy output integrated over flare 



Table 8: Ratio of synchrotron to SSC or EC flare amplitude expected when a physical 
parameter varies. 



Parameter varied 


SSC/Synch 


EC/Synch 


b 


fa 1 


> 1 




> 1 


^ 1 


B 


fa 1 


< 1 



-28- 



w 

'e 

o 

0) 



X 

3 



E 
o 



X 

3 



w 

M 

'e 

o 

0) 



X 

3 



T ' 

Monitor 



t.n I n f !i t II I uL Li 




1996 



1998 



2000 



2002 



2004 



2006 



2008 




1996.94 



1996.98 



1997.02 
YEAR 



1997.06 



1997.1 



Fig. 1. — X-ray (2-10 keV) data on different time-scales. In tfie upper panel, the arrows 
show the times of superluminal ejections and the line segments perpendicular to the arrows 
show the uncertainties in the times of ejection. 
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Fig. 2. — Optical (R-band) data on different time-scales. 
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Fig. 3. — Radio (14.5 GHz) data on different time-scales. 
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3C 279 

VLBA/7mm 




-1 -2 -3 



Fig. 4. — Sequence of VLBA images at 7 mm during 2001. The images are convolved with 
the beam of the size 0.38x0.14 mas at PA = —9°. The global peak over all maps is 17.24 
Jy/Beam. The contour levels are 0.25, 0.354, 0.5, 0.707, 90.51 % of the global peak. The 
angular scale given at the bottom is in milliarcseconds (mas). 
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Fig. 5. — Sequence of VLBA images at 7 mm during 2002. For details see caption of Figure 
4. 
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Fig. 6. — Sequence of VLBA images at 7 mm during 2003. For details see caption of Figure 
4. 
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Fig. 7. — Sequence of VLBA images at 7 



mm 



during 2004. For details see caption of Figure 
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Fig. 8. — Sequence of VLBA images at 7 mm during 2005. For details see caption of Figure 
4. 
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Fig. 9. — Sequence of VLBA images at 7 mm during 2006. For details see caption of Figure 
4. 
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3C 279' • 




1996 1998 2000 2002 2004 2006 

Epoch 

Fig. 10. — Angular separation from the core vs. epoch of all knots brighter than 100 mJy 
within 2.0 mas of the core. The black lines indicate the motion of each knot (denoted by a 
given color of data points) listed in Table 2. A knot is identified through continuity of the 
trajectory from one epoch to the next. The diameter of each symbol is proportional to the 
logarithm of the flux density of the knot, as determined by model fitting of the VLB A data. 
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Fig. 11. — Variation of X-ray flux, optical flux, radio flux and position angle of the jet from 
1996 to 2008. The circled data points in the bottom panel are the epochs shown in Fig. [T^l 
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3C 279 




2 1 0-1-2-3 -4 



Fig. 12. — VLBA images at one epoch during each year of 11-year monitoring. The images 
are convolved with the beam of the size 0.38 x0.14 mas at PA = —9°. The map peak is 
17.0 Jy/Beam. The contour levels are 0.15, 0.3, 0.6, ...,76.8 % of the peak. The angular 
scale given at the bottom is in milliarcseconds (mas). The circled points in Fig. [11] (bottom 
panel) correspond to these images. 
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YEAR 

Fig. 13. — Light curve of the VLBA core region at 43 GHz. The jagged hne through the 
data points is drawn solely to aid the eye to follow the variability. Statistical and systematic 
uncertainty in each measurement is difficult to determine accurately, but is typically 10-20%. 
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log [Frequency (Hz)] 

Fig. 14. — Result of application of the PSRESP method to the light curves. PSD of the 
observed data at high, medium and low frequency range are given by the solid, dashed 
and dotted jagged lines respectively while the underlying power-law model is given by the 
dotted straight line. Points with error bars (open squares, solid circles and asterisks for 
high, medium and low frequency range respectively) correspond to the mean value of the 
PSD simulated from the underlying power- law model (see text). The errorbars are the 
standard deviation of the distribution of simulated PSDs. The broadband PSD in all three 
wavelengths can be described by a simple power law. 
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Fig. 15. — Success fraction vs. slope for all three PSDs. The success fractions indicate the 
goodness of fit obtained from the PSRESP method (see text) . 
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Fig. 16. — Discrete cross-correlation function (DCCF) of the optical, X-ray, and radio mon- 
itor data. The time delay is defined as positive if the variations at the higher frequency 
waveband lag those at the lower frequency. 
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Fig. 17. — Variation of X-ray/optical time lag across overlapping 2 intervals from the begin- 
ning of the first year to the end of the second, except for the bottom right panel, for which 
the interval is indicated more precisely. Notice the big change between 1998-99 and 2000- 
2001, when the time delay went from X-ray leading optical to the opposite sense. There is a 
major change between 2001-02 and 2003-04 as well, when the time delay went from optical 
leading X-ray to the opposite sense. In the last 4 years (bottom right panel) the correlation 
became weaker but maintained the negative time delay. 
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Fig. 18. — Cross-correlation function of the X-ray light curve and the position angle of the 
jet. Changes in the position angle lead those in the X-ray flux by 80 ± 150 days. 
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Fig. 19. — Smoothed X-ray and optical light curves. Curves correspond to summed flux 
after modeling the light curve as a superposition of many individual flares. Thick horizontal 
strips in the X-ray light curve in 2000 correspond to epochs when no data are available. 
Flare pairs hsted in Table 4 are marked with the respective ID numbers. 
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Fig. 20. — Cross-correlation of the X-ray and 43 GHz core light curves. Changes in the 
X-ray flux lead those in the radio core by I30I45 days- 
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Fig. 21. — Simulated synchrotron (solid curve) and SSC (dashed curve) flares. Here all 
flares are created by an exponential rise and decay in the magnetic field B [B, Nq, and R 
are functions of distance along the jet). In the bottom panels, •jmax is increased linearly with 
time causing the SSC flux to peak ahead of the synchrotron flux. Flare amphtudes have 
been scaled such that they can be seen on the same plot. Normalization in the two upper 
panels is the same and that in the two lower panels is the same. 
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Fig. 22. — Simulated synchrotron (solid curve) and SSC (dashed curve) flares. Here all flares 
are created by an exponential rise and decay in the magnetic field Nq [B, Nq, and R are 
functions of distance along the jet). In the bottom panels, ■jmax is increased linearly with 
time causing the SSC flux to peak ahead of the synchrotron flux. Flare amphtudes have 
been scaled such that they can be seen on the same plot. Normalization in the two upper 
panels is the same and that in the two lower panels is the same. 
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Fig. 23. — Segments of real light curves (optical: asterisks, X-ray: open circles) over ~ 
100 day intervals similar to the length of the simulated light curves. Optical flux density is 
multiplied by R band central frequency (4.7 x 10^^ Hz) to obtain the flux. The simulated 
light curves in Figures EH and [22] have profiles similar to flares in the actual light curves. 



